Mesoscopic fluctuations of the Coulomb drag 
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We consider mesoscopic fluctuations of the Coulomb drag coefficient po in the system of two 
separated two-dimensional electron gases. It is shown that at low temperatures sample to sample 
fluctuations of po exceed its ensemble average. It means that in such a regime the sign of pr> is 
random and the temperature dependence almost saturates po ~ 1/VT. 

PACS numbers: 73.23.-b, 73.50.-h, 73.6I.-r 
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When two electronic layers are brought close to each 
other to form a bi-layer system, a current flowing through 
one of the layers(the active layer) is known to induce 
a voltage Vd — PdI in the other (passive) layer [Q-Q. 
The effect, which is called the drag, was first predicted 
theoretically |l],|^] in the model where the carriers in 
the two spatially separated layers interacted via long- 
range Coulomb interaction. Experimentally the Coulomb 
drag was first observed in a three-dimensional electron 
gas layer while the current was driven through a two- 
dimensional electron gas (2DEG) H). Subsequent exper- 
iments studied the effect in 2DEG bi-layers jf|, electron- 
hole ||, and normal-metal-superconductor systems ||. 
More recently, the drag was studied in the 2DEG bi-layer 
system in high magnetic field . 

The quantity, which is studied theoretically is the 
transconductance ap>. To the lowest non- vanishing or- 
der in the interlayer interaction it is proportional to the 
drag coefficient pp> {<Ji is the Drude conductance of the 
i-th layer) 
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As a function of temperature the observed ap, roughly 
follows the quadratic law ap ~ T 2 , although the ratio 
ctd/T 2 deviates from the constant value 0,0. 

The T 2 dependence of the Coulomb drag coefficient 
follows from the Fermi liquid phase space argument. To 
create a current in the passive layer, it is necessary to 
create a pair of electron-like (filled states with energy 
greater than the Fermi energy e > ep) and hole- like 
excitations (empty states e < ep) in a state with non- 
zero momentum. The energy and momentum of the pair 
come from an electron in the active layer, which is mov- 
ing with the driving current. In each layer, the scatter- 
ing states are limited to the energies of order T relative 
to the Fermi level, which gives two powers of T to the 
drag coefficient. However, the momentum is transferred 
equally to electrons and holes, therefore in the case of 
electron-hole symmetry the drag of the electrons cancels 
that of the holes. Thus the effect is non-zero only due to 
the electron-hole asymmetry. Similarly, the asymmetry 
is necessary for the electron and hole system in the active 



layer to have non-zero total momentum in the first place. 
The asymmetry can be expressed as a derivative of the 
density of states v and/or the diffusion constant D with 
respect to the chemical potential p. This can be obtained 
rigorously in the diagrammatic formalism p0| . For the 
case of diffusive layers the disorder-averaged transcon- 
ductance is 



[°D 



(hT) 2 



d 



h 3 g 2 {nd) 2 \dp 



[vD] 



In 



2T' 



(1) 



where for simplicity we take the layers to be identical, 
so that they have the same chemical potential, diffu- 
sion constant and the dimensionless conductance g = 
25.8kCl/Rn- The logarithm is cut at the scale To = Dn/d 
and k = 2ire 2 v is the inverse Thomas- Fermi screening 
length. 

Such effects of the electron-hole asymmetry are well- 
known, for instance the thermopower in disordered elec- 
tronic systems jll]] or adiabatic pumping |l2). As these 
effects are due to the electron-hole asymmetry, the av- 
erage quantities are small, since each derivative with re- 
spect to the chemical potential brings one power of the 
Fermi energy Ep to the denominator. On the other hand, 
the typical energy scale of mesoscopic effects is the Thou- 
less energy Ep — T1D/L 2 (L is the sample size), which 
is much smaller than the Fermi energy. Therefore the 
effects mentioned above exhibit mesoscopic fluctuations, 
much larger than the average. 

In this Letter we show that the mesoscopic fluctuations 
of the Coulomb drag coefficient can indeed be larger than 
the average Eq. ([!]), even if the electron systems in both 
layers are good metals (j> 1). To characterize the fluc- 
tuations, we calculate the average square of the (random) 
transconductance. The result of lengthy albeit straight- 
forward calculations shown below is given by 
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where the numerical factor 7 = 1.0086 is the value of the 
integral 
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' = H / dx / 1 f C ^ Xl ^ 2 [J (xi)J (x 2 ) + J2 (11)^2(^2)] 
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where Ji(x) and Ki(x) are the Bessel functions fl3|| . Here 
t v <C (Et) -1 is the dephasing time. This result is valid in 
the most relevant regime L v = */ Dt v < L and nd ^ 1. 
If K<i < 1, then the average square of the conductance 
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i D l g 4 with the coefficient of order unity. 

In what follows we first discuss the experimental conse- 
quences of our results, then explain it qualitatively and 
finally give the rigorous calculation. The fluctuations 
Eq. (J2J) depend on temperature only through the dephas- 
ing time 1/t v w T/g and at low enough temperatures 
they should dominate the behavior of the transconduc- 
tance. Therefore the T 2 decrease of <td should at some 
small temperature T* be almost saturated at a sample- 
dependent value. Let us estimate T* for the samples used 
in existing experiments using the reported parame- 
ters of the samples. Collecting the numerical factors, we 
write the ratio of the square of the average transconduc- 
tance Eq. ([!]) and the averaged square Eq. (|J) as 
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We take the interlayer spacing to be d = 200 A [[l||7]] ; the 
screening length in GaAs is k _1 = 100A; the Thouless 
energy is given by Et = g/(2-KvL 2 ); and the dephasing 
Tlng/g Then we estimate T* as the 
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temperature at which the ratio (0\d) 2 /(°z)) is equal to 
unity 
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where the Fermi energy in the samples ||,|7|] Ep ~ 60A', 
the electron density n = 1.5 x 10 n cm~ 2 , the size of the 
sample L ~ 200/im, and the conductance is calculated 
from the sheet resistance of the sample R = 1051/n. 

The estimated T» is lower than the temperature range 
for the existing data jijj^], therefore there is no trace of 
the fluctuations Eq. (|) in the data. However, if one takes 
a dirtier sample, with the sheet resistance, for instance, 
lkil/O, then the estimate for T* becomes 2K and the 
effect of the fluctuations becomes observable. To push 
T* even higher, the sample size can be also reduced. 

Let us now explain Eq. (|^) qualitatively. First, con- 
sider the lowest temperature regime T <§; Et, so that the 
sample is effectively zero-dimensional (0D). The meso- 
scopic fluctuations of the usual conductance are universal 
5a ~ e 2 /Ti. The transconductance is associated with the 
interlayer interactions, thus possessing additional small- 
ness. The value of such smallness can be estimated by 
the Golden rule argument which is comprised by (i) phase 



volume; (ii) matrix elements; (iii) electron-hole asymme- 
try (dependence of the density of states on the energy) . 
Matrix elements in 0D do not depend on energy and 
give smallness 1/g 2 JT^]. Therefore the phase volume is 
limited by temperature only, which gives the factor T . 
Finally, the electron-hole asymmetry d^ihw^)) is the 
random quantity with the typical value 1/Et- Putting 
everything together, we arrive to the estimate 
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FIG. 1. A. The diagram for the transconductance in 
the lowest non-vanishing order in the interlayer interac- 
tion. The solid lines are the exact Green's functions of 
the non-interacting electron system (in the presence of dis- 
order) and the wavy lines are the disorder averaged inter- 
layer Coulomb interaction propagators Eq. (Q). B. The RPA 
scheme for the calculation of the interaction propagators. The 
dashed lines denote the disorder. The indices i, k, I = 1, 2 in- 
dicate the layer. 

At higher temperatures the averaging should be per- 
formed by dividing the sample into patches of the size 
L v x Lip, since on larger scales the phase coherence is 
destroyed. The contribution of each patch (^^(L^)) is 
approximately the same and they can be simply com- 
boncd as a network of random resistors 
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Now, (5a™ s (L v )) can be found similar to Eq. (|3|) with 
two important differences: (i) the fluctuations of the den- 
sity of states are summed from the scale of order T, 
rather than Et{L v ). This suppresses the fluctuations 
in each layer by the factor yj Et(L v )/T. (ii) The matrix 
elements become energy dependent on the energy scale 
larger than Et, decreasing with the transmitted energy 
uj as \M\ 2 ~ 1/ui 2 . As a result, the transmitted energy is 
limited by u> ~ Et(L v ) — l/r v , rather then by T, so in 
the estimate of the phase volume we should replace T 2 
by T/t v . So, we find 
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Finally, to estimate the total magnitude of the transcon- 
ductance fluctuations we substitute Eq. @ into Eq. (|) 
to obtain 



This estimate yields the same result as Eq. (g|) up to 
numerical factors and the dependence on n,d. 

Our results suggest the following picture of the 
Coulomb drag. If one starts measuring the drag co- 
efficient at high T and proceeds by lowering the tem- 
perature, then at first the transrcsistance will decrease 
roughly as T 2 , as follows from Eq. ([j]). At the tempera- 
ture T„, estimated above, the transresistance will appear 
to saturate (or> oc 1/y/T), as the fluctuations Eq. (g) will 
start to dominate. The particular value of the prefactor 
will be sample dependent and, what is more important, 
will have random sign. If the temperature will be de- 
creased further, then at very low temperatures T < Et 
the sample will effectively become zero dimensional and 
the T 2 decrease will be restored (also with a random co- 
efficient), so that at T = the drag coefficient vanishes. 

Let us now present the calculation. The electrons in 
both layers interact via the Coulomb interaction. The in- 
teraction propagators corresponding to the dynamically 
screened Coulomb interaction can be obtained within the 
RPA scheme (see Fig. IB) using the Green's functions 
of non-interacting electrons. For our purposes we only 
need the propagator of the interlayer interactions, which 
is given by (here we set the layers to be identical, so 
that they have the same density of states v and diffusion 
coefficient D) 
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The transresistance in the disordered two-layer system 
can be expressedin terms of the exact Green's functions of 
non-interacting, disordered electron system and the inter- 
action propagators Eq. (||). To the lowest non- vanishing 
order in the interlayer interaction the transresistance is 
given by the diagram Fig. 1A. The left and right triangles 
correspond to the two layers in the system and the wavy 
line is the interlayer interaction propagator Eq. (|^). As 
the electron Greens's functions now depend on disorder, 
this ob is a random quantity and its moments should be 
averaged over disorder. Before averaging, the expression 
for <7d corresponding to the diagram on Fig. ^ can be 
written as 

D 4V J 2ir \duj 2T J 12 23 34 41 ' V ; 

where the indices indicate the spatial coordinates. Points 
1, 2 belong to one layer and 3, 4 to the other. The trian- 
gular vertices T a are given by 
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The exact electronic Green's functions used in Eq. (g) 
can be written in terms of the exact wavefunctions of the 
system as 
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where j labels the exact eigenstates of the system and 6j 
are the exact eigenvalues. 

The known result for the averaged transconductance 
Eq. (@) is obtained @ by averaging the triangular ver- 
tices T a independently for each layer (for the effect of the 
correlated disorder see Ref. |l6|]). However, as we dis- 
cussed above, in the intermediate temperature range the 
fluctuations exceed the average, and the temperature de- 
pendence saturates. To characterize the fluctuations we 
average the square of the transconductance. The meso- 
scopic fluctuations of the interaction propagators can be 
neglected, because they produce only the small fluctuat- 
ing coefficient in Eq. (|7|). Therefore, we have to average 
the product of two triangular vertices T a (in the same 
layer; note that in this temperature regime Eq. (jScl) does 
not contribute). The corresponding diagrams are shown 
in Fig. 2. After averaging, each diagram contains six dif- 
fusons and one Hikami box (see Fig. 3). The calculation 
of the arising 14-dimensional integral is greatly simplified 
by the following observation. The dominant contribution 
to the frequency integrals Eq. (|8a| ) comes from the region, 
where the external frequencies [which are the frequencies 
of the interaction propagators Eq. (^)] uii + ui2 — nd/r v 
and io\ — 102 — 1/ T ^ so that the frequency difference is 
small compared to the sum. The momentum integral is 
dominated by the region \/t v < Q < nd/r^. Since we are 
in the intermediate temperature regime T > Et, the en- 
ergy transfer u> is smaller than the temperature, and the 
vertices Eq. (||) can be expanded in the inverse temper- 
ature. Now the dimensional analysis gives the resulting 
expression for the fluctuations Eq. (j^), which depends 
on temperature only through t v as we discussed above. 
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The numerical coefficient can now be obtained by per- 
forming the integration without further approximations. 
The factor 7 comes from the angle integration and the 
numerical factor in Eq. ( pbj ) from the integration over the 
small frequency difference. 
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experiments |3Hj. For the samples used in |J,|7[ we esti- 
mated the crossover temperature T» « 0.2K, which was 
below the temperature range used in these experiments. 
To observe the effect of the fluctuations, one needs to 
take a dirtier sample of smaller size. Then T* can be 
equal to several Kelvin, and the saturation of o\d to a 
value with the random sign can be observed. 

Finally, we notice that the Coulomb drag coefficient 
may also be presented as the product of two random 
numbers po — &iet2, where aj., 0,2 characterize each layer. 
If the disorder is corelated, the average (0.1112) appears, 
which leads to the results of Ref. |l6| . In this respect, 
the results of Ref. |l6| are just a particular manifestation 
of the mesoscopic fluctuations of po, discussed in this 
Letter. 
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FIG. 2. Difhison contribution to the average of the product 
of two triangular vertices. Cooperon contribution is obtained 
by interchanging of vertices p and j in one of the triangles 
(the direction of arrows should be changed respectively) 






FIG. 3. Hikami box. 

In conclusion, we have described the mesoscopic fluc- 
tuations of the Coulomb drag coefficient or the transcon- 
ductance. The fluctuations are characterized by the aver- 
age square of the random, disorder dependent transcon- 
ductance Eq. (0). Compared to the averaged transcon- 
ductance Eq. §\§ the fluctuations Eq. (||) are determined 
by the Thouless energy, rather than the Fermi energy, as 
the average. Therefore, there exists an intermediate tem- 
perature regime, where the fluctuations are greater than 
the average, which results in the weak (1/yT) tempera- 
ture dependence of the transconductance in this regime. 
Moreover, in this regime the <jd is a random, sample 
dependent quantity, so that the sign of the measured 
value is also random. Since the average transconductance 
Eq. (jl]) grows as T 2 , at higher temperatures (T > T*) the 
fluctuations are small and the the measured ub is roughly 
equal to the average. This was the case in the existing 
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